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ABSTRACT 

We show that a fast-cooling synchrotron spectrum has a photon index harder 
than -1.5, as long as one introduces the evolution of magnetic field in the emis- 
sion region during the course of emission, an important physical process hitherto 
ignored in theoretical modeling. By introducing a reasonable decay law of mag- 
netic fields, we model the prompt emission spectrum of GRBs, and found that 
the low-energy photon index below the injection frequency is typically around -1 
in the time scale of seconds, well consistent with the data. The model photon 
spectrum can be fit with the well-known "Band" function. All these suggest that 
fast-cooling synchrotron radiation is the dominant mechanism to power GRBs. 



Subject headings: gamma-ray burst: general 
thermal 



radiation mechanisms: non- 



The radiation mechanism of gamma-ray bursts (GRBs), the most luminous explosions 
in the universe, remains un-identified after 45 years since their discovery in late 1960s. Syn- 
chrotron radiation o f electrons accelerated in relativistic shocks has be en suggested as the 
leading mechanism (IMeszaros et al.l Il994l ; iDaigne fc Mochkovitchl Il998l ). Strong magnetic 
fields in the emission region lead to rapid cooling of electrons (fast cooling). Traditionally it is 
believed that the fast cooling electrons give ri se to a photon sp ectrum with a power law spec- 
tral index -1.5 below the injection frequenc y (pari et al. 119981). while obseryationally the typ 



ical low-energy photon spectral index is -1 (IBand et al. 



2011; 



Nava et al. 



199 



3; 



Preece et al.ll2000t IZhang et al. 



2011). As a result, fast cooling synchrotron mechanism has been disfavored 
200o( ), and slow cooling due to decay of magnetic field from the shock front 



(Ghiselli ni et all._ -, , - - . 

( Pe'er k, Zhang! 2006 ) or slow heating due to turbulent acceleration ( Asano k, Terasawa 2009 
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Murase et al.l 120121 ) have been invoked in order to inte rpret the GRB spectrum with syn- 



chrotron radiation. Other mechanisms were proposed (IRees &: Meszarod 120051 ; iPe'er et al. 



20061 ; iBeloborodovl |2010| ; lLazzati et al.l l2013t Kumar fc Panaitescul 120081 ) , but observations 
and theoretical argume nts suggest that they may not be the dominant mecha nism to power 
the majority of GRBs Jzhang et aliboij IVurm et al. l l2013l : IPiran et al.ll2009h . 



A typical GRB promp t emission spectru m is a smoothly connected broken power law 
called the "Band function" ( Band et al.lll993l ). The high-energy photon spectral index (with 
a typical value /3 ~ —2.2) can be easily interpreted as emission of fast-cooling electrons with 
an injected power law distribution Q(7e) oc 7~ p , given a typical injected electron spectral 
index p ~ 2.4. The low-energy photon spectral index (with typical value a ~ —1), on 
the other hand, has defied all known t heoretical models. The photosphere model typically 
predicts a = +0.4 jBeloborodovll2010h. In order t o reproduce the data, some special jet 



structure has to be introduced (jLundman et al.ll2013f ). For synchrotron radiation, an ironclad 
value a = —1.5 (jSari et al.lll998l ) has been widely believed in the community. 



This well-known index a = —1.5 can be derived from a simple argument. Let us consider 
a continuity equation of electrons in energy space (d / dt) (dN e / dj e ) + {d / dj e ) [ r j e (dN e /d n f e )] = 
Q(7e, t), where dNg/dje is the instantaneous electron spectrum of the system at the epoch 
t, and Q(j e ,t) is the source function above a minimum injection Lorentz factor / ~f m of the 
electrons. For synchrotron radiation, the electron energy loss rate is 



7e 



67rm e c 



oc 



(1) 



where ax, Tn e , and c are Thomson cross section, electron mass, and speed of light, respec- 
tively, and B is the strength of magnetic fields in the emission region. For fast cooling, 
electrons are cooled rapidly to a particular energy j c (t) (cooling energy) below the injection 
energy j m . In the regime 7 C < 7e < 7m, one has Q(j e ,t) = 0. Also consider a steady state 
system (d/dt = 0), then one immediately gets dN e /d^ e oc 7~ 2 , i.e. the electron spectral 
index is p = 2. The specific intensity of synchrotron spec trum would have a spectra l index 
s = (p — l)/2 = 1/2 (with the convention F v oc v~ s ) ( iRybicki fc Lightmanl Il979l ). The 



photon spectral index (defined as dN^/dE^ oc E", where E 7 is the photon energy, and iV 7 is 
the photon number flux) would then be a = 



7 : 

1 + s) 



■1.5. 



The above argument relies on a crucial assumption of a steady state, which is achieved 
when a constant B is invoked. In a GRB problem, since the ejecta launched from the central 
engine is expanding in space, the magnetic field strength cannot be preserved as a constant. 
In the rest frame of a conical jet, flux conservation indicates that the radial (poloidal) 
magnetic field component decreases as Bl oc r~ 2 , while th e transverse (toroidal) magnetic 
field component decreases as B' t oc r -1 (jSpruit et al.l 120011 ). As a result, at a large radius 
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from the central engine where 7-rays are radiated, one has a toroidal-dominated magnetic 
field with B' oc r _1 . Here r is the distance from the central engine. Considering other effects 
(e.g. magnetic dissipation), the decay law may be described by a more general form 



, -b 

r 



B ' {r) = B '°[~o) • (2) 

We investigate a generic problem of electron fast cooling in a decaying magnetic field de- 
lineated by Eqj2j and study the synchrotron emission spectrum. This model does not specify 
a particular energy dissipation mechanism or particle acceleration mechanism, a nd hence, can 



apply to a variety of GRB prompt emission models such as internal shocks (IDaigne et al. 



2011 ) and internal collision-induced magnetic reconnection and turbulence (jZhang Sz Yan 
201ll ). We consider a toy box that contains electrons and a co-moving magnetic field B', 
which moves relativistically towards the observer with a bulk Lorentz factor T. The rela- 
tivists electrons are accelerated into a power-law distribution Q(y e , t') = Qo(t')[y e / r y rn (t')]~ p 
(for 7 e > 7 m (i')) of a slope p and continuously injected into the box at the injection rate 
Rinj(t') = f°° Q{pj e ,t')drf e , where t' is the time measured in the co-moving fluid frame. Here 
Rm]{t') St' gives the number of electrons injected into the box during the time interval t' and 
f + St'. 

Electrons undergo both radiative and adiabatic cooling. In the rest frame that is co- 
moving with the relativistic ejecta , the evolution of the Lorentz factor y e of an electron can 



be described by (jUhm et al.l 120121 ) 



<1 ( 1 A o"t B ,2 _ 1 ( J_\ dlnn e ^ 



dt' \7e/ 6nm e c 3 \7 e / dt 

For a conically expanding toy box, we take the comoving electron number density n e oc r~ 2 , 
which gives dhxn e = — 2c?lnr. We divide the injection function Q(j e ,t') into small divisions 
in time space t' and also in the energy space j e , and numerically follow cooling of each group 
of electrons (between [f, t' + St'] and [y e , 7 e + Sj e ]) individually using Eqj3l We then find 
the instantaneous global electron spectrum d~N e jd^ e of the system at the current epoch. 

We consider four models with different decay indices b in Eqj2j Model [a] takes a constant 
co-moving magnetic field B' — B' Q — 30 G. This implies that there should be no change in the 
volume of the box. Thus for this model we drop out the adiabatic cooling term from Eqj3j 
As shown in Column 1 of Fig.l and Fig.2, this model gives the familiar electron spectrum 
dNe/d^e cx 7~ 2 below 7 m , and the well-known photon spectrum F u oc v~ x l 2 in the fast cooling 
regime. Model [b] takes ro = 10 15 cm and b = 1.0 in EqfSJ This is the case of free conical 
expansion with flux conservation (no significant magnetic dissipation). As shown in Column 
2 of Fig.l and Fig.2, the electron spectrum and the photon spectrum both harden with 
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time. At 1.0 s after injecting the first group of electrons, the global electron energy spectral 
index deviates significantly from the p = 2 nominal value, and hardens to around p ~ 1. 
The corresponding photon spectrum is nearly flat (F u oc u°), which correspond s to a photon 



index ~ —1, the typica l low-energy photon index observed in most GRBs (IPreece et al. 



20001 ; IZhang et al.ll20!ll ). This gives a very simple interpretation to the mysterious low- 
energy photon index, which has challenged theorists for decades. In Columns 3 and 4 of 
Fig.l and Fig.2, we present Models [c] and [d], for which a steeper decay index b = 1.2 
and b = 1.5 are adopted, respectively. These cases may mimic the c ases when significan t 



magnetic dissipation occurs during the course of synchrotron radiation (IZhang &: Yanll201ll ). 



As shown in Fig.2, both models also give spectra that are consistent with the observations. 

In order to understand the physical origin of such an effect, in Fig. 3 we decompose the 
^obs = 1-0 s instantaneous electron spectrum into the contributions of 10 injection time slices, 
each lasting for 0.1 s. For the constant B' case (Fig.3a), one can see that as the electrons 
age, they tend to distribute more narrowly in logarithmic energy space, so that the electron 
number per energy bin increases. This is because in the fast cooling regime, as time elapses, 
the original electrons with a wide range of energy distribution tend to cool down to a narrow 
range of cooling energy defined by the ages of the electrons in the group, which are very close 
to each other at late epochs. Above j m , the electron energy distribution remains unchanged 
with time, since it is always determined by the same injection rate and cooling rate. 

The B' decay cases show a more complicated behavior. The distribution of each group 
of electrons still shrinks as the group ages. However, since at early epochs the magnetic 
field was stronger, it had a stronger cooling effect so that for a same injection time duration 
(0.1 s), initially it had a wider spread in energy at a given age (which can be noticed by 
comparing the 0.1 s electron spectrum for Model [a] and [b] in Fig.l). The later injected 
electrons are cooled in a weaker B' field, so that their initial spread is narrower. After the 
same shrinking effect due to cooling pile- up, the groups injected in earlier time slices have 
a wider electron distribution than the constant B' case. Also the electron spectrum above 
the injection energy, even though possessing a same spectral index, has a normalization 
increasing with time due to progressively less cooling in a progressively weaker magnetic 
field. These complicated effects all work in the direction to harden the spectral index, as 
seen in Figures 3b, 3c and 3d. For a steeper .B'-decay index (e.g. b = 1.2 and b = 1.5), the 
late-time injection occurs in an even weaker magnetic field, so that slow cooling is possible. 
This results in accumulation of electrons around the minimum injection energy 7 m , so that 
a sharper break in the electron energy distribution is achieved. 

The model predicts that the low-energy spectrum below the injection frequency v m is 
curved, due to the complicated cooling effect as delineated in Fig. 3. Most GRB detectors 
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have a narrow band pass so that below the peak energy (typically a few hundred keV), there 
are at most 2 decades in energy. Nonetheless, in the detector band pass, the observed spectra 
are usually fit by a Band function, with the low energy spectral index a ~ — 1. Usually tim e 
resolved spectral analyses are carried out with a time bin in seconds ( IZhang et al.ll201ll ). 
This is the typic al time scale of the slow variability component in most GRB light curves 
( iGao et al.ll2012l ). We therefore focus on the 1 s and 3 s model spectra. We truncate these 
spectra in a narrow band (5 keV - 5 MeV) and compare them with the empirical Band 
function fits (Fig. 4). One can observe that most of our model spectra are consistent with 
the Band function with the correct low-energy spectral indices. 

Outside the band pass, our model predicts an asymptotic value of the low-energy electron 
energy spectral index of p a = (66 — A) /(6b — 1), which is 2/5 for 6 = 1. This is seen in the 
numerical results of the models (lower panels in Fig.l), and can be derived analytically (see 
Appendix). According to the simple relationship s = (p — l)/2, one gets s a = —3/ (126 — 2), 
which is -0.3 for 6 = 1 (or F v oc u 3 ). In reality, due to the contribution of the 1/3 segment 
of the individual electron spectrum, which becomes significant when p approaches 1/3 from 
above, the asymptotic photon spectrum limit is softened. In this case, it is about -0.2. This 
corresponds to a photon index of -0.8, which is much harder than the nominal value -1.5. 
The spectrum is somewhat softer in a narrower band 5 keV - 5 MeV, but is close to -1 
(Fig.4). 

Our model also predicts that during the very early epochs of a pulse (e.g. the first 0.1 s, 
0.01 s, etc), the spectrum should be softer. However, usually there are not enough photons 
during these early epochs to give a statistically significant spectral analysis. Unusually bright 
GRBs with distinct pulse structures may be used to test this prediction in the future. 

The models presented here are valid when inverse Compton cooling is not important, 
which requires a highly magnetized e.jecta, as has been suppo rted by recent observations 
(IZhang fc Pe'erl l2009t IZhang et al.ll2011t lYonetoku et al.ll2012l ). The example models pre- 
sented here only invoke a constant injection function. GRB lightcurves show an erratic 
behavior, so in reality Q and 7 m should be time-dependent. This would give rise to richer 
spectral features. 
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A. Asymptotic Low Energy Shape of Electron Spectrum 



The asymptotic low-energy spectral index can be derived analytically for Eq.3. For 
prompt GRB emission, it is reasonable to adopt a constant Lorentz factor T, so that one has 
r = ct'T. We first solve a simpler equation by dropping out the adiabatic term, i.e. 



where 



d ( 1 

Gnm P c 



at 



,-2b 



Bl (cr/r ) 



-2b 



We then find the solution of electron Lorentz factor at any time t'- (> £•) 



7e(0 = 



a 



,, 1-26 _ ,,1-26 

1 - 2b I j 



(Al) 
(A2) 

(A3) 



where 7 e (^) is the electron Lorentz factor at an initial time t\. For b > 1/2, i' ^> t^, and 



leit'j) <C 7e(0, this solution gives 7e (£'-) oc ^ 2b 1 . We then get 



*7e(**) « if - 2 ^ oc [ 7e (<J.)] ^ ^. 



(A4) 



For a constant injection rate i?i n j, we have 0~iV e oc 8t[. Thus, we have an asymptotic behavior 
of the global electron spectrum as follows. 

5N, 



OC 



26-2 
' 26-1 



(A5) 



Now we consider the full Eq.3 that includes the adiabatic term. For B(r) = B (r/r ) b 
and r = ct'T, Equation 3 can be written as 



d f 1 

This equation has an analytic solution 

7e(0 = t'~ 2/3 



,,-26 , 2 / 1 

at + 



3t> V7 



3a 



1-66 



/(l-66)/3 



(A6) 



(A7) 



where C is the integration constant of the differential equation, to be determined by the 
initial condition; 7e (t-) at time t\. The electron's Lorentz factor 7e (^) at a later time t'- is 
found to be 

,,-2/3 o 

h, , f ,(l-66)/3 /(l-66)/3\ / A8 x 



7«(0 = t'- 2/3 



For b > 1/6, t'j > ^ and tf 3 ^) < tf /3 7e(t-), this solution gives 7e (t$) oc t'f 2/3 tf b ' 1)/S . 
A variation in 7 e (^) results only from St^ for the instantaneous (i.e., at a fixed time ^) global 
electron spectrum; <? Te (t$) oc t^ 2/3 ^ {6fe " 4)/3 5^. This gives the asymptotic behavior of the 
global electron spectrum 

^yoc^W^.)]-^, (A9) 

where we have again assumed a constant injection rate i?i n j. Therefore we have the asymp- 
totic low-energy electron spectral index p a = (66 — 4)/ (66 — 1). 
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Fig. 1 . — The co- moving frame fast cooling electron energy spectrum evolution as a function 
of time t'. The injected electrons have a power law distribution Q(j e ,t') = Qo(t')('j e / r y m )~ p 
above a minimum injection Lorentz factor 7 m = 10 5 , with a power law index p = 2.3 for 
7e > 7m- We take a co-moving magnetic field B' = B' (r/r )~ b with B' Q = 30 G and r = 10 15 
cm. Four models are investigated: [a] 6 = 0; [b] 6 = 1; [c] 6 = 1.2; and [d] b = 1.5. For all the 
models, a constant injection rate = J°° Q(7 e ,t')^7e = 10 47 s _1 is adopted, with both Q 
and 7 m as constants. The electron injection into the box begins at r = 10 14 cm. The ejecta 
is assumed to be moving towards the observer with a Lorentz factor F = 300, and the burst 
is assumed at a cosmological redshift z — 1. For each model, the instantaneous electron 
spectra at four different epochs since the beginning of electron acceleration are calculated. 
The four epochs in the observer's frame are: 0.1 s (black), 0.3 s (blue), 1.0 s (red), and 
3.0 s (green). For each epoch, the sharp cutoff at low energies corresponds to the "cooling 
energy" of the system, which is defined by the strength of the magnetic fields and the age 
of the electrons. Given the same B' field at r = 10 15 cm, the B' field is stronger at earlier 
epochs, so that electrons undergo more significant cooling initially. One can see that the 
cooling energy is systematically lower than that of the constant B' case (Model [a]) as b 
steepens (Models [b], [c] and [d]). The lower panel of each model shows the local electron 
spectral index as a function of electron energy j e . For Model [a] (constant magnetic field), 
the electron spectrum shows the well-known broken power law, with the spectral indices p+1 
and 2 above and below the injection energy, respectively. For other models, even though the 
index above 7 m remains the same, the index below 7 m is much harder. At later epochs (e.g. 
3 s spectra), the index approaches an asymptotic value p a = (66 — 4)/ (66 — 1). 



- 10 - 




Fig. 2. — The synchrotron emission spectra of electrons with energy distribution p resented 
in Fig.l. The full synchrotron spectrum of each electron ( iRybicki fc Lightmanlll979h is taken 
into account. The observed spectra are calculated by considering the Lorentz blueshift and 
cosmo logical redshift. While the constant B' case (Model [a]) gives rise to the familiar 
F u oc z/ -1 / 2 spectrum, the decaying B' cases (Models [b], [c] and [d]) all give rise to much 
harder spectrum below the injection break v m . For the spectra in the seconds time scale (1 
s - red; 3s- green), the low-energy spectral index is nearly flat, consistent with the typical 
observed photon index -1. Lower panels show local spectral slopes as a function of observed 
frequency. 
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Fig. 3. — Decomposition of electron spectrum at 1 s in the observer's time. In order to see 
the contributions of electrons injected at different epochs, the electrons are grouped into 10 
slices in injection time, each with a duration of 0.1 s. The contributions of each group to 
the instantaneous electron spectrum at 1 s are marked in different colors. Older groups are 
cooled down further towards lower energies, so from left to right, curves with different colors 
denote the electron energy distribution of the electron groups injected from progressively 
later epochs, with a 0.1 s time step. Dashed curves are the summed total of all electrons. 
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Fig. 4. — A comparison of our 1 s and 3 s model spectra (solid) with the empirical Band 
function fits (dashed) for all four models in a narrower band pass from 5 keV to 5 MeV. It 
is seen that the model spectra can mimic the Band function spectra well. The plotted Band 
function parameters are the following: Model [a]: a = —1.5, (3 = —2.1, E = 1800 keV for 
both 1 s and 3 s; Model [b]: a = -1.2, (3 = -1.95, E = 480 keV for 1 s, and a = -1.15, 
(3 = -1.95, E = 220 keV for 3 s; Model [c]: a = -1.15, f3 = -1.92, E = 390 keV for 1 s, 
and a = -1.13, (3 = -1.92, E = 190 keV for 3 s; Model [d]: a = -1.1, (3 = -1.9, E = 330 
keV for 1 s, and a = -1.05, (3 = -1.8, E = 100 keV for 3 s. 



